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ABSTRACT 

The  radiated  acoustic  field  from  a  point-driven, 
fluid-loaded  plate  was  determined  several  years  ago. 
However,  current  operational  problems  force  a  re- 
examination in  the  light  of  certain   acoustic  anomalies. 

The  types  of  waves  that  propagate  in  a  solid  plate 
are  reviewed  for  their  relative  importance. 
Experimental  models  were  constructed  in  an  anechoic 
acoustic  tank  in  an  attempt  to  model  the  aforementioned 
anomalies . 

It  is  found  that  the  effect  of  a  plexiglass  plate 
placed  at  the  boundary  of  a  method-of- images  geometry  is 
that  of  a  phase  delay  only.  Therefore,  mode  conversion 
from  compress ional  waves  to  flexural  waves  is  thought  to 
be  insignificant. 

A  second  experiment  is  described  in  which  a  plexi- 
glass plate  was  driven  with  a  small  longitudinal 
vibrator.  Experimental  radiation  patterns  suggest  that 
of  a  quasi-piston.  Plexiglass  is  an  acoustically  slow 
media,  therefore  the  coincidence  effect  is  not  possible. 
Further,  its  absorption  of  acoustic  waves  make  it  a  poor 
choice  for  experiments  where  wave  propagation  is 
important . 
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I.  INTRODUCTION 

A.   THESIS  INTRODUCTION 

The  experiments  on  which  this  thesis  is  based  were 
attempts  to  physically  model  a  sound  source  mounted  on 
the  ice  in  the  Arctic  using  a  sheet  of  plexiglass  and 
piezoelectric  transducers  in  an  anechoic  water  tank. 
Specifically,  the  sound  field  radiated  into  the  fluid 
was  of  interest. 

There  are  many  wave  types  that  propagate  in  a  solid 
and  along  a  solid-fluid  interface.  Each  was  examined 
for  its  relative  importance,  so  that  attention  could  be 
concentrated  on  those  that  are  more  important  to  the 
radiation  of  sound  into  the  water  from  a  sound  source 
mounted  on  the  ice.  The  Ewing  and  Press  [Ref.  1]  method 
of  solution  covers  wave  propatation  in  ice  far  from  the 
source.  Stein  [Ref.  2]  considered  radiation  from  sound 
sources  mounted  on  the  ice. 

Two  experiments  were  performed.  The  first  was  an 
investigation  on  the  effect  of  a  plate  placed  at  the 
free  boundary  of  the  classic  me thod-of- images  (Lloyd's 
mirror)  geometry.  Waves  emanating  from  a  source  near  a 
pressure-release  boundary  will  develop  an  interference 
pattern  with  waves  reflected  from  the  boundary  that  can 
be  detected   by  a  receiver.    It  is   hypothesized  that  a 
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plate  placed   at  the  boundary  will   affect   the  inter- 
ference pattern. 

In  the  second  experiment  performed,  the  source  was 
attached  to  the  ice  in  an  attempt  to  duplicate  the 
"acoustic  mode,"  named  by  Stein  [Ref.  2,  p. 100],  of  the 
Arctic  ice  sound  propagation.  The  theoretical  develop- 
ment of  the  chapter  uses  the  Mindlin-Timoshenko  thick 
plate  equations  to  explain  the  effects  of  coincidence 
and  fluid  loading  on  the  radiated  sound  field  from  a 
point-driven,  fluid-loaded  plate. 

B .   BACKGROUND 

The  classical  two-dimensional  theory  of  flexural 
waves  in  elastic  bars  in  vacuo,  using  Lagrange's 
equation,  was  solved  in  1889  by  Lord  Rayleigh  under  the 
assumption  that  the  wavelengths  were  long  with  respect 
to  the  thickness  of  the  plate.  Thus,  the  equations  are 
invalid  for  high  frequency  waves  or  for  sharp 
transients.  There  are  two  effects  that  were  neglected: 
rotatory  inertia  and  transverse  shear  deformation. 
Later,  however,  Rayleigh  introduced  into  the  equations 
the  effect  of  rotatory  inertia,  but  because  transverse 
shear  deformation  was  not  included,  this  correction  was 
too  large.  Timoshenko  [Ref.  3]  included  the  effect  of 
transverse  shear  deformation  in  a  study  of  flexural 
waves  of  bars. 
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Mindlin  [Ref.  4]  in  a  study  of  flexural  waves  in 
plates  included  the  effects  of  both  rotatory  inertia  and 
shear  deformation  in  his  now  classic  paper  in  which  are 
developed  the  so-called  Mindlin-Timoshenko  thick  plate 
equations.  Dyer  [Ref.  5],  using  the  Mindlin-Timoshenko 
equations,  showed  the  effects  of  transverse  shear  to  be 
important  even  for  thin  plates.  Skudrzyk  et  al  [Ref*  6] 
developed  a  theory  of  flexural  waves  in  a  driven  plate 
similar  to  that  of  the  theory  of  reverberation  in  room 
acoustics . 

In  1966,  David  Feit  [Ref.  7]  calculated  the  acoustic 
pressure  radiated  into  a  fluid  by  an  infinite  plate 
driven  by  a  point  force,  using  the  Mindlin-Timoshenko 
thick  plate  equations.  His  solution  used  the  saddle 
point  method  of  steepest  descent  to  evaluate  the  pres- 
sure integral.  Stuart  [Ref.  8]  explored  the  solution  of 
the  Mindlin-Timoshenko  equations  with  then- new  insights 
concerning  the  so-called  leaky  waves  emanating  from  the 
plate.  He  asserted  that  previous  discrepancies  in  the 
saddle  point  technique  solution  (a  coincidence  frequency 
lower  than  predicted,  and  small  but  unaccounted  for 
lobes  in  the  far-field  pattern  below  coincidence)  were 
due  to  neglected  poles  occurring  near  the  saddle  point, 
called  leaky  wave  poles,  and  that  these  influence  the 
far-field  pressure. 
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Creighton  [Ref.  9]  contradicted  these  findings, 
demonstrating  that  "...  the  far-field  contributions 
associated  with  .  .  .  the  leaky  waves  may  be  altered 
merely  by  a  change  in  the  integration  path  ..."  and 
that  the  contributions  decayed  exponentially  with  depth 
into  the  water  and  did  not  represent  physically  meaning- 
ful free  modes,  except  for  subsonic  surface  waves. 

As  the  theory  of  waves  in  a  plate  developed,  so  did 
the  theory  of  wave  propagation  in  the  ice.  The  theory 
of  propagation  of  elastic  waves  in  a  floating  ice  sheet 
was  given  by  Press  and  Ewing  [Ref.  1],  by  considering 
the  ice  as  an  infinite  plate  and  using  classical  wave 
theory.  Press  and  Ewing  [Ref.  10]  then  extended  their 
results  to  include  the  effects  of  the  air.  A  further 
discussion  of  their  work  as  a  whole  is  found  in  Ewing, 
Jardetzky,  and  Press  [Ref.  11]. 

In  a  recent  study  Stein  [Ref.  2]  used  the  Ewing  and 
Press  method  of  solution,  but  solved  the  problem  with 
the  aid  of  modern  computer  algorithms.  The  bulk  of  the 
study  was  that  of  the  propagation  speeds  of  the  various 
modes.  His  data  taken  from  Arctic  and  freshwater  lake 
ice  experiments  corraborate  with  the  predicted  transient 
arrival  times  and  propagation  losses . 

C.   WAVES  IN  A  SOLID 

There  are  several  wave  types  that  may  propagate  in  a 
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solid  and  others  that,  occur  at  a  solid-fluid  interface. 
A  study  of  the  applicability  of  each  type  found  in  the 
literature  was  done  and  is  included  here  to  give  the 
reader  a  sense  of  wave  types  that  are  relevant  to  the 
problem  of  modeling  a  sound  source  mounted  on  the  ice 
generating  a  sound  field  in  the  underlying  water  by  the 
excitation  of  the  ice. 

Propagation  of  waves  in  solids  can  be  distinguished 
from  that  in  fluids  by  the  solid's  ability  to  restore 
shear.   This  leads  to  two  wave  equations: 


S*       3 


"at* 


(l.la) 


for  the   scalar  displacement   potential  which   describes 
longitudinal  motion,  and 


v2  ;--A 


32 


3- 


(1.1b) 


for  the  vector  displacement  potential  which  describes 
shear  wave  motion.  The  potentials  <j>(x,y,z,t)  and 
t(x,y,z,t)   are   related  to  the  particle  velocity  vector 

.  A  A  A 

u  =  uxi  +  uvj  +  uzk   by 


u 


V<fr  +  Vx<ji 


(1.2) 
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An  excellent  depiction  of  some  of  the  different 
types  of  acoustic  waves,  along  with  their  speeds,  is 
given  by  Ver  and  Holmer  [Ref.  12,  p.  273]  and  it  is  from 
this  source  that  Figure  1.1  is  taken. 

Longitudinal  waves  are  waves  that  propagate  along 
the  direction  of  compression  and  dilatation  in  the 
plate.  The  thickening  and  thinning  of  the  plate  seen  in 
Fig. 1.1a  is  a  consequence  of  a  non-zero  Poisson's  ratio. 

Shear  waves  are  solutions  to  the  additional  wave 
equation  that  governs  waves  in  solids : 

It 


t«*-^,  m 


cs 


where    c«2  ■  —  as  in  Fig.  1.1b 
s     P 


ij>  ■  vector  potential  of  particle  displacement 
in  a  solid. 


Two  types  of  shear  waves  are  found  in  the 
literature:  SH  and  SV  waves  are  the  designations  given 
them  in  texts  on  seismological  wave  propagation.  SV, 
shear  waves  polarized  vertically,  denotes  waves  for 
which  the  shearing  motion  is  perpendicular  to  the  plane 
of  the  plate.  It  is  this  wave  that  is  depicted  in  Fig. 
lb.  SH,  shear  waves  horizontally  polarized,  denotes 
waves  for  which  the  shearing  motion  is  parallel  to  the 
plane  of  the  plate. 
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Figure  1.1. 


Wave  Types  in  Solids  and  Their  Speeds 
(After  Ver  and  Holmer) 


Y  «  Young's  modulus,   which   relates   the   stress   to  the 

strain  produced, 
a   ■  Poissons  ratio, which  relates  lateral  expansion  to  its 

longitudinal  shortening   due  to  a  compressive  stress, 
G  *  the  shear  modulus,  which  relates   shearing  stress  to 

shearing  strain  produced, 
I  =»  the  moment  of  inertia, 
K  =  a   torsional  stiffness  factor,   which  relates  a  twist 

to  the  shearing  strain  produced, 
p  =  the  solid's  density, 
t  -  the  solid's  thickness,  and 


u> 


the  radian  frequency  of  the  motion 
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Torsional  waves,  Fig.  lc,  result  from  a  twisting 
force  applied  to  a  solid.  This  wave  can  be  generated  in 
bars,  where  the  y  and  z  dimensions  are  similar. 
Torsional  waves  are  not  be  easily  excited  in  a  plate; 
they  will  not  be  further  considered  in  this  study. 

Flexural  waves,  or  bending  waves,  are  a  more  complex 
motion  than  the  previous  ones  considered,  because  of  the 
complicated  stresses  set  up  when  a  plate  is  bent.  The 
bending  of  a  plate  compresses  the  plate  on  the  inside  of 
the  bend  and  expands  the  plate  on  the.  outside.  Due  to 
this  compression  and  expansion,  a  Poisson's  ratio  effect 
is  introduced. 

As  with  any  wave,  the  energy  of  a  flexural  wave  is 
composed  of  two  parts,  potential  energy  and  kinetic 
energy.  The  potential  energy  is  made  up  of  two 
deflections:  flexural  deformation,  and  shear 
deformation.  The  kinetic  energy  also  has  two 
contributions,  linear  (or  trans latory)  inertia,  and 
rotatory  inertia. 

The  classical  two-dimensional  theory  of  flexural 
motions  in  thin  elastic  plates  leads  to  the  wave 
equation  [Ref.  4,  p.  31] 


DVSjr  +  ph  |^f  =0  (1.3) 

Yh3 
where    D  ■  .  2  j  * — 2T    *    an<^  here 
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h  =  plate  thickness 

w  =  the  plate  particle  motion 
This  equation  relates  only  the  flexural  deformation  to 
the  linear  (or  trans latory)  inertia  as  the  two  opposing 
forces.  The  Mindlin-Timoshenko  equation  relates  the 
linear  and  rotatory  inertia  as  the  driving  forces  to  the 
flexural  and  transverse  shear  deformation  as  the  restor- 
ing forces.  (This  was  what  was  referred  to  earlier  as 
"correcting"  for  rotatory  inertia  and  transverse  shear.) 
The  Mindlin-Timoshenko  equation  is 

(1.4) 

(*2  "  jr  &)  (dv2  -  £i  fa)     w  t  phtf^  w  =  o 

where 


w2     Y 


12   TTT+a7) 
h   =  plate  thickness 


o   32 
The  -!rr  . .  ■>  term  is  the  contribution  of  transverse  shear 


Dh  3  %2 

The  ^ 37T  term  is  the  contribution  of  rotatory 

inertia. 


Note  that  if  these  terms  are  neglected,  the  classical 
equation  (Eq.  1.3)  results. 
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The  Mindlin-Timoshenko  theory  is  used  when  analyz- 
ing "thick"  plates,  i.e.,  plates  whose  thickness  exceeds 
Xs/20  ,  where  Xs  is  the  shear  wave  length.   [Ref.  13,  p. 

215] 

Of  key  importance   to  this  study   is  the   effective 
phase  speed  for  flexural  motions  as  a  function  of 
frequency.    Ross,  assuming   a  harmonic  disturbance  of  a 
homogeneous  plate   in  vacuo,   gives  an   equation  for  the 
flexural  wave  speed   c^   as 


(1_  r  (JL.) )*(<*) *  +  (rp+i)(£-)2(^)2  -  (£-)2  =  o   (1.5) 
*   "p    cp  t>    p     rp 


where 


Tp  =  the  plate  shear  parameter 

=  2.65  (1  +  |  a   +  |  a2) 

flp  =  the  reference  angular  frequency, 

.  cp   12 


a)  ■  the  angular  frequency  of  the  flexural  wave 

Cp  ■  the  plate  longitudinal  wave  speed,  given  in 
Fig.  1.1a. 

Equation  1.5  can  be  solved  by   applying  the   quad- 


ratic equation  to  solve  for  the  variable   (-°-k)   .   It  is 

CP 

a  difficult  function  of  the  other  parameters,  but  can  be 
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represented  accurately  in  terms  of  approximate  solutions 
at  low,  intermediate  and  high  frequencies. 

At  low   frequencies,   the   flexural   wave   speed   is 
approximated  by 


cbl 


i 


0)2Yt2 
12p(l-az) 


(1.6) 


which  is  that  given  in  Fig.  l.ld,  and  represents  the 
flexural  wave  speed  in  thin  elastic  plates,  given  by  the 
classical  wave  equation,  Equation  1.3.  This  equation  is 
valid  when 


rP<%'    K<1 


In  the   intermediate  frequency  range,   the   flexural 
wave  speed  is  approximated  by 


cbi 


=  c 


bl 


d+'i:'2) 


1/2      •. 

-  2  TL 


-1  1/2 


1  "  'p  V2 


(1.7) 


In  the  high  frequency  range,  the  flexural  wave  phase 
velocity  predicted  by  the  Mindlin-Timoshenko  theory 
tends  to  the  Rayleigh  wave  velocity  [Ref.  13,  p.  248], 
which  is  a  value  slightly  less  than  the  shear  wave 
velocity.     The    expression   for   the   high   frequency 
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flexural  wave  speed  is 


cbh  =   0.76(1  +  |  a)    c 


s 


(1.8) 


This  equation  is  valid  when 

Heretofore  we  have  discussed  waves  propagating  in  a 
solid.  Waves  can  also  propagate  along  the  boundary  of  a 
solid-fluid  interface  in  much  the  same  manner  as  do  sur- 
face waves  on  the  ocean.  These  waves  are  Rayleigh  waves 
(for  a  half-solid),  Lamb  waves  (for  a  plate),  Stonely 
waves,  Love  waves,  and  Crary  waves. 

Rayleigh  waves  " .  .  .represent  elastic  perturbations 
propagating  near  the  free  boundary  of  a  solid  and  decay- 
ing with  depth"  [Viktorov,  Ref.  14,  p.  1] .  Figure  1.2 
is  a  depiction  of  the  instantaneous  deformation  associ- 
ated with  a  Rayleigh  wave. 


—  —  — .  — —  i 1 

—  _  _«__.—  —  —'  —  '—'  —  —>  —  —  —  _  _  — . 


Figure  1.2   Depiction  of  a  Rayleigh  Wave  (after 
Viktorov) 
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Rayleigh  waves  propagate  along  the  surface  of  a 
solid  at  a  free  boundary,  such  as  an  air-ice  boundary  or 
an  air-plexiglass  boundary.  The  wave  decays  exponen- 
tially in  the  z-direction  into  the  solid,  and  if  the 
decay  constant  is  not  large,  the  motion  of  the  wave  may 
appear  at  the  lower  (in  this  case,  water)  boundary. 
Furthermore,  the  Rayleigh  wave  is  supersonic  with 
respect  to  the  fluid,  therefore  this  motion  would  radi- 
ate sound  into  fluid  (see  Chapter  III  for  an  explanation 
on  the  importance  of  relative  sound  speeds  to  radiation 
into  the  fluid.) 

The  discussion  of  Rayleigh  waves  is  included 
because,  as  mentioned  above,  in  the  high  frequency  limit 
flexural  waves  degenerate  to  Rayleigh  surface  waves  for 
a  plate  in  vacuo.  However,  Stein  shows  that  at  low  fre- 
quencies there  is  no  coupling  of  the  Rayleigh  wave  with 
the  flexural  wave  there  would  be  for  a  plate  in  a 
vacuo1.  Moreover,  the  types  of  sources  being  considered 
would  not  generate  Rayleigh  waves,  hence  they  will  not 
be  further  considered. 

Lamb  waves  "...  refer  to  elastic  perturbations 
propagating   in   a   solid   plate  (or   layer)   with   free 


^here  is  a  growing  debate  among  the  cognoscenti 
concerning  the  importance  of  the  Rayleigh  wave  in  Arctic 
ice  [Ref.  15]. 
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boundries,  for  which  displacements  occur  both  in  the 
direction  of  wave  propagation  and  perpendicular  to  the 
plane  of  the  plate."  [Ref.  13,  p.  67].  They  are  also 
called  the  "normal  modes"  of  a  plate*  and  can  be  sym- 
metrical (8)  or  antisymmetrical  (a)  relative  to  the 
center  plane  of  the  plate.  Note  that  the  zero-order 
symmetric  and  antisymmetric  modes  are  the  longitudinal 
and  flexural  waves,  respectively.  Figure  1.3  depicts 
the  zero-order  Lamb  waves.  These  two  modes  are  of  key 
importance  to  this  study. 


W 


Figure  1.3   Lamb  Wave  Propagation  (after  Viktorov) 

Another  wave  that  can  exist  at  an  interface  between 
two  media  is  a  Stonely  wave.  According  to  Ewing, 
Jardetzky  and  Press  [Ref.  11,  p.  112]  they  are  modified 
Rayleigh  waves  that  exist  when  the  shear  wave  velocities 
in  the  two  media  are  nearly  equal.  Stein  [Ref.  2.  p. 
69]  states  that  they  can  also  exist  along  the  boundary 
between   an    elastic  medium   and   a   high    density, 
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acoustically  slow  fluid  or  elastic  medium,  and  that  the 
Stonely  wave  is  the  high  frequency  asymptote  of  the 
first  ant i- symmetric  mode  along  the  ice-water  interface. 
That  is,  they  are  Rayleigh  waves  along  the  solid-fluid 
interface.  However,  absorption  in  Artie  ice  is  high  at 
the  frequencies  for  which  Stonely  waves  exist  [Ref.  2, 
p.  190]  and  therefore  the  Stonely  wave  is  not  important 
to  this  study. 

Two  other  wave  types  concerning  propagation  of  sound 
in  ice  are  described  in  the  literature.  Love  waves  and 
Crary  waves.  Both  are  shear  waves.  Love  waves  are  SH 
waves  that  are  generated  due  to  discontinuities  in  the 
ice.  Stein  [Ref.  2  p.  100]  notes  that  "...  unless 
there  is  a  discontinuity  ....  they  will  not  radiate 
nor  can  they  couple  to  other  modes."  Discontinuities  in 
the  ice  are  not  considered  in  this  study,  therefore  Love 
waves  will  not  be  considered. 

Crary  waves  are  a  type  of  SV  wave  discovered  by 
Crary  in  seismic  experiments  on  a  floating  ice  island. 
They  were  generated  after  an  ice  "event",  a  small  ex- 
plosive charge  mounted  in  the  ice.  According  to  Ewing, 
Jardetzky  and  Press  [Ref.  11,  p.  299],  this  wave  has  the 
following  characteristics: 

"1.   phase  velocity  is  near  the  speed  of  compression- 
al  waves  in  the  plate; 
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2 .  travel  time  is  intermediate  between  that  for  P 
(compressional)  waves  and  that  for  SH  waves; 

3.  principal  [amplitude]  is  [observed]  on  longitud- 
tudinal  horizontal  seismograph;  amplitudes  are 
much  smaller  on  vertical  and  practically  absent 
on  transverse  horizontal; 

4.  frequency  is  almost  constant,  increasing  very 
slightly  with  time." 

Stein  [Ref.  2,  p.  93]  again  discusses  these  waves, 
but  has  found  the  attenuation  to  be  at  least  0.1  dB/m 
in  the  horizontal  for  the  Arctic.  Thus,  these  waves 
will  not  be  observed  at  appreciable  distances.  Further, 
because  of  their  small  vertical  amplitude  and  velocity 
less  than  that  of  the  fluid,  (characteristic  3  of  above) 
they  do  not  impart  acoustic  energy  to  the  underlying 
fluid.  For  these  reasons,  Crary  waves  will  not  be 
further  considered. 

So  far  in  this  chapter,  the  types  of  waves  that  pro- 
pagate in  a  solid,  namely  Arctic  ice,  have  been  dis- 
cussed without  regard  to  their  origin.  That  is,  they 
have  been  assumed  to  be  plane  waves,  far  from  their 
source.  As  these  waves  propagate,  they  may  be  coupled 
into  the  surrounding  medium,  creating  acoustic 
radiation. 

Additional  radiation  may  be  generated  when  the  waves 
are  excited  by  a   localized  force.    This  is  a  result  of 

27 


the  flexural  motion  having  a  spatial  gradient  of  dis- 
placement in  the  ice,  which  plane  waves  do  not.  Stein 
fRef .  2  p.  963  nas  named  this  node  of  acoustic  radiation 
the  "acoustic  mode . H 

Figure  1.4  is  a  depiction  of  the  first  antisymmetric 
mode,  the  flexural  wave,  in  a  solid.  At  this  point,  far 
from  the  source,  it  can  be  seen  why  the  flexural  wave 
does  not  radiate2:  it  travels  subsonically  with  respect 
to  the  water  and  there  is  cancellation  of  radiation  from 
adjacent  positive  and  negative  pressure  peaks.  (In 
actuality,  because  of  geometric  spreading,  there  is  in- 
complete cancellation  and  some  radiation  does  exist.) 
However,  when  these  waves  are  excited  by  a  localized 
force   (Fig.  1.5)   there  is  incomplete   cancellation  and 


Figure  1.4   Flexural  Wave  Propagation 


2A  flexural  wave  below  coincidence  will  not  radiate. 
The  coincidence  frequency  is  that  frequency  above  which 
the  plate  vibrations  can  couple  to  the  fluid.  See 
Chapter  3  for  a  further  discussion. 
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radiation  occurs.  Note  that  this  acoustic  mode  is  still 
a  flexural  wave.  A  further  discussion  of  this  mode  found 
in  Chapter  III. 


Vacuum 


/|\ 


Fluid 


Figure  1.5   Flexural  Wave  Propagation  From  a  Point 
Source  (after  Stein) 

D .   SUMMARY 

In  this  chapter,  the  many  types  of  waves  found  in  a 
solid  (ice)  have  been  discussed.  As  was  indicated 
earlier,  only  a  few  have  real  importance  in  the  Artie, 
and  therefore  modeling  in  the  lab  should  be  attempted 
only  for  these.  Stein  [Ref.  2  p.  107]  has  concluded 
that  the  three  modes  of  propagation  of  importance  are 
the  first  symmetric  and  antisymmetric  modes,  and  the 
acoustic  mode.  All  other  aforementioned  waves  are 
neglected  because  either  they  are  highly  attenuated  in 
the  Arctic  at  the  frequencies  of  concern  or  there  is  a 
lack  of  a  source  with  which  to  generate  them. 
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II.   A  LLOYD'S  MIRROR  EXPERIMENT 

A.  INTRODUCTION 

In  this  chapter,  an  experiment  will  be  discussed 
that  tested  the  following  hypothesis:  a  plate  intro- 
duced at  an  air-water  interface  will  influence  the 
acoustic  field  in  the  water  half-space.  More  simply, 
the  reflection  of  a  spherical  wave  from  an  ideal 
boundary  develops  an  interference  pattern  with  the 
original  incident  wave.  It  is  hypothesized  that  the 
plate  will  effect  this  Lloyd's  mirror  pattern. 

B.  THE  METHOD  OF  IMAGES 

The  method  of  images  is  a  well-known  technique  for 
predicting  the  radiation  pattern  of  an  acoustic  source 
near  a  fluid  interface.  For  continuity,  a  brief  de- 
scription of  the  theory  of  the  acoustic  Lloyd's  mirror 
will  be  presented.  The  following  description  is  based 
on  that  given  by  Sanders  and  Coppens  [Ref.  16]. 

If  a  source  of  spherical  waves  is  placed  on  the  z- 
axis  a  distance  d  below  a  pressure-release  boundary,  as 
shown  in  Fig.  2.1,  the  acoustic  radiation  from  this 
source  in  the  absence  of  the  boundary   can  be   specified 

as        -   =  *    eJUt-kr_)  (2.oi) 


*    r. 


where 


■_  =   J(h-d)2  +  r2 
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As  a  result  of  the  boundary,   there  must  be  an  addi- 
tional wave  present  to  describe  the  reflection: 


p_r  =  -  £-   eJUt-kr+)  (2.02) 


where  r+  =  J(h+d) 2  +  R2 
The  total  field  is,  then:  (2.03) 


£s  +  £r  =  £ 


^r_  r+       J 


The  acoustic  field  will  show  strong  interference  effects 
because  the  two  waves  p_s  and  p_r  must  travel  along 
different  paths  to  reach  the  receiver. 

The  geometry  is  simplified  if  the  assumption 

r>>dcos6  (2.04) 

is  made.   Then,  r+  and  r_  can  be  replaced  by 

r+  -  r+Ar  (2.05) 

r_  «  r-Ar 
Ar  »  dsin8 
with  2.04,  2.03  becomes 

E(r,e.t)  ^e^^D^.e^j 

r        r 


31 


Figure  2 . 1   Lloyd ' s  Mirror  Geometry 

Thus  when  a  simple  harmonic  acoustic  wave  in  the  the 
water  is  incident  on  a  plane,  perfectly  reflecting  pres- 
sure release  boundary,  the  reflected  wave  is  a  mirror 
image  of  the  incident  wave  which  leads  to  an  interfer- 
ence pattern  in  the  medium  below  the  boundary. 

This  interference  pattern  may  be  observed  by  dif- 
ferent methods.  The  receiver  can  be  placed  at  a  fixed 
horizontal  range  from  the  source  and  moved  downward  to 
show  the  nodes  and  antinodes  from  the  pattern.  It  might 
also  be  placed  at  a  fixed  slant  range  from  the  source 
and  swept  through  an  arc  about  the  source.   Or,  it  might 


32 


be  held  at  a  constant  depth  and  moved  away  from  the 
source  in  horizontal  range.  Lastly,  and  want  was  done 
in  the  experiment,  the  source  and  receiver  can  be  held 
in  a  fixed  position  and  the  frequency  of  the  source 
swept,  observing  the  response  of  the  receiver.  The 
interference  of  the  incident  and  reflected  waves  produce 
maxima  and  minima  in  the  frequency  response  of  the 
receiver. 

C.   THE  EFFECT  OF  THE  PLATE 

The  effect  of  the  plate3  on  the  incident  wave  can  be 
best  seen  by  reverting  to  simple  ray  theory.  Consider 
the  following  problem:  an  incident  wave  in  water,  re- 
presented by  ray  p_i»  impinges   on  a  plate   of  different 

characteristic  impedance  r   ,  as  shown  in  Fig .  2.2.   The 

2 

ray  (now  p_a)  is  refracted  as  it  traverses  into  the 
plate.  It  reflects  from  the  top  of  the  plate  at  z  =  H 
(now  gfc) ,  traverses  the  plate  and  then  is  refracted  back 
into  the  water  as  pr  .  Note  that  while  the  wave  is  in 
the  plate,  it  is  considered  as  a  compressional  (longitu- 
dinal) wave . 


3An  infinite  plate  forced  by  an  incident  sound  wave 
at  any  frequency,  even  below  its  coincidence  frequency, 
will  radiate  sound,  although  a  bending  wave  of  the  plate 
itself,  below  coincidence,  will  not.  [Ref.  12,  p.  280] 
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Figure  2.2   Ray  Traversing  Through  a  Plate 


The  four  waves  may  be  represented   as  (considering   only 
two  dimensions): 

£i  =*  Zi  eJUt-*  sineix-k  cos  8^)  (2.1) 

Da  =»  P.  •J<»t"*  sin8ax-k  cos8az)  (2.2) 

i.a    — a         2         2 

£b  "  Zb  eJUt-k  sin8bx+k  cos8bz)  (2.3) 

p  ■  p   eJ(wt-k  sin8rx+k  cos8rz)  (2.4) 

*-r   — '        1         1 


The  boundary  conditions  are: 

at    z  »  0,  continuity  of  pressure  requires 


£i(x,o)  +£r(x,o)  -£a(x,o)  +  ^(xjO) 


(2.5) 
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Also  at  z=0,  the  condition  of  continuity  of  normal 
velocity,  requires  (taking  into  account  direction  of 
propagation) 

Ui(x,o)cos9r  -  Uj-Cx/OjcosSj- 

(2.6) 
=   ua(x,o)cos9a-ub(x,o)cos9b 

The  boundary  at  z  ■  H  is  pressure  release,  therefore 

£a(x,H)  +  p^UfH)  =  0  (2.7) 

For  each  plane  wave,  we  have 

£  =  ±r  (2.8) 

u 

With  Snell's  Law,     sin  9^   sin  9*   ,  (2.9) 

cl        c2 

and  since  the  angle  of  incidence  equals  the  angle  of 
reflection,  we  have 

9j_  ■  6r  =  9X   and   9a  =  9b  =  92  (2.10) 

If  each  value  of  u  in  2.6  is  replaced  with  the 
appropriate  value  of  p/r  ,  and  from  2.10,  subscription 
are  changed,  then: 

SfafiifeiU.oH^U.o)]  -52£i2  [£a(x,o)-Eb(x,o)](2-12) 

*1  r2 


or 


£i(x,o)-£r(x,o)  -£l|2|£2  (Ea(x#o)-£b(xfo))   (2.13) 

i  o  COS  \j   1 
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Apply  2.1  through  2.4  to  2.5  ,  and  drop  the  e3ut  , 


-jk  sine  x  -jk  sine  x 

(£i+ir>e    l  l      ■  (la+£b>e    2     2         (2.14) 


and  do  the  same  to  2.13, 


(p  p^)e-jklsin8ix  .  r^  cose^  (p-p^-  jk28ine2x(2-15 , 

^■"  r   cos  v 


Apply  the  boundary  condition  at  z=H  (Equation  2.7),   and 


-jk  sine  x 
dividing  through  by  e  2     2   gives 


pae-Jk2cose2H  +  ^e+jkzeosezH  =  Q  (2>16) 

i 


Algebraic  manipulation  of  2.14  -  2.16  then   leads  to  the 
reflection  coefficient  R  : 


p     (i  _  £l  cos92-)  _  n  +  rj,  cose 2)  e-j2k2cose2H 

,  ^£  _  *■ r?  cose  j  J    * r ?  cose  t  j  

£i    Ci    £i  cos9?  ^  (-.         r\    cose^  -  j2k2cose2H 

*>     r  cose  '  ^     r   cose  ' 

2      1  2  1  (2.17) 


where  e   and  e   are  related  by  Snell's  Law  (.9),   and 
1       2 


cose2   =   Jl    -    (^)2    sin^  (2.18) 


■>Mt 
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D.   VALIDATION  OF  EQUATION  2.17 

In  the  limit  as  H+O  ,  Equation  2.17  reduces  to 
R  =-l  .  This  is  to  be  expected,  because  when  H+0  the  r 
layer  disappears  and  the  incident  wave  is  2  reflected 
off  a  purely  pressure-release  boundary,  where  the  re- 
flected wave  is  180°  out  of  phase  with  the  incident,  and 
hence  R  =  -1  . 

In  the  limit  as  9  +0  ,   that  is,   normal   incidence, 

l 

Equation  2.17  reduces  to: 

„  -rt  cosk2H  +  jr2sink2H  (2.19) 

— *■    rj  cosk2H  +  jr2sink2H 

Comparing  Equation  2.19  with  Equation  6.12  of 
Kinsler  et  al  [Ref.  17] 

R  _  [1~  ft]cosh2H  +  Jtrl"  f^)sinh2H  (220) 

(1+  f^COSh2H  +  *(ff  +  ft)sinh2H 

where   r   is  characteristic  impedance  of  the  fluid  above 
3 

the  place.    In  the   limit  r  +0  ,   that  is,   assuming   a 

3 

vacuum  above  the  plate,  2.20  reduces  to  2.19. 

Finally,  comparison  of  2.12  with  Brekhovskikh  [Ref. 
18]  yields  additional  validation.   His  coordinate  system 
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which  reduces  to  [Ref.  18,  Equation  5.18] 

R  m   y„  ±   V1?  ez^2d 

-    1  +  V   V    z:ia2d 
2  3    12e 

where   V12  =  f^f*  .   V23  =  §*=§*  (2.25) 


Assuming   Z   =0   (in  vacuo),   then  V   =  -1 ,   and   .27 
1  12 

reduces  to 

tt   <*--fy  -  u^)e2j^d  (2.26) 

-   (1  +  §*)  -  (1  -  |l)e2ja2d 

2  2 

changing  to  previous  notation  and  noting  that  d  =  -H  , 
equation  2.26  becomes  2.17: 


t.         r_L  cos  8  o  \  /-.    r_j_  cos  9  ?  >>   -zjk2cos92H 

I    r?  cose  i  ^  "  ^ r?  cose  t  j 

—    fi  +  £_L  cos8?  ^  r-i    £j.  cos_8_2>,i   -zjk2cos82H 

*>■    r  cose  *  ^    r  cose  ' 

2       1  2  1 


E.   VALIDATION  OF  THE  PLANE  WAVE  ASSUMPTION 

Thus  far,  the  reflection  coefficient  of  an  oblique 
incident  plane  wave  off  a  finite  layer  has  been  calcu- 
lated and  validated  against  previously  known  theory  and 
limiting  cases.   However,  one  question  still  needs  to  be 
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addressed  before  this  equation  is  put  to  use.  That  is, 
the  validity  of  the  plane  wave  assumption  made  above, 
Equation  2.8,  needs  to  be  established. 

The  solution  for  a  point  source  near  a  boundary  is 
addressed  in  two  works.  Officer's  [Ref.  19]  solution  is 
in  the  form  of  three  integrals  for  each  of  the  three 
regions.  However,  the  integrals  are  difficult  to 
compute.  Officer  does  give  "...  the  condition  under 
which  the  plane  wave  reflection  coefficient  may  be  taken 
as  a  valid  approximation  for  the  reflection  of  an  acci- 
dent spherical  wave.   This  condition  is  that 


7-» <<  1   "  (2.27) 

C,^     ...o.i  3/2 


«  Ri  [^-4  -  sin^e] 

2 


where  R^  is  the  distance  from  point  of  reflection  to  the 
receiver. 

The  second  work  that  addresses  this  problem  is 
Brekhovskikh.  From  Brekhovskikh,  the  geometric  optics 
approximation  (in  which  a  spherical  wave  is  relected  off 
a  layer  with  the  same  reflection  coefficiient  as  a  plane 
wave)  is  valid  when  the  distances  of  the  source  and  the 
receiver  from  the  boundary  are  large  compared  with  the 
wavelength.  [Ref.  18  p.  256]  In  the  experiment  run, 
this  assumption  is  valid  above  f  =  50kHz. 
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Lastly,  for  completeness,  the  extent  of  the  near 
field  should  be  calculated.  From  Ziomek  [Ref.  20],  the 
maximum  radial  extent  r  of  the  near  field  of  the  source 
is  given  by: 


r  > 


JLjl  (2.28) 


where 

R   is   the  maximum  radial   dimension   of  the  transmit 

transducer. 

For  the  experiment  done  here,  the  following  values 
were  used  in  equations  2.27  and  2.28: 

ci   =  1481  m/s 

c2  ■  1992  m/s 

6X  m   42.27° 

R^  =  0.5  m 
Equation  2.27  requires  that  f  >>  14.8  kHz  for  the  plane 
wave  assumption  to  be  valid.  From  2.28,  for  the  radia- 
tion from  the  transmitting  LC-10  to  be  in  the  far-field 
prior  to  reaching  the  plate  requires  f  <  1.01  MHz. 
Therefore  for  the  frequency  range  of  interest  (50-100 
kHz)  both  conditions  are  satisfied. 

F.   APPLICATION  OF  EQUATION  2.17 

Upon  close  inspection  of  Equation  2.17,  it  can  be 
seen  that  the  reflection  coefficient  is  a  phasor  of  unit 
magnitude  with   a  phase  delay   that  is   a   function   of 
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frequency.  This  is  intuitively  true  because  of  the  con- 
ditions imposed  on  the  problem:  the  boundary  condition 
at  the  top  of  the  plate,  Z  =  H  ,  is  a  pressure  release, 
perfectly  reflecting  plane;  the  medium  is  assumed  to  be 
lossless.  Since  it  is  perfectly  reflecting,  there  is  no 
refracted  wave  at  the  upper  boundary  and  there  is  no 
loss  of  signal.  If  the  medium  is  lossless  there  is  no 
loss  of  signal  as  the  wave  travels  through  the  plate. 
Thus,  the  reflection  is  of  unit  magnitude  with  a  phase 
angle,  that  is  R  =  l*e^T. 

The  phase  angle  t  of  R  is,  from  Equation  2.17, 


t  =  tan 


-,  (1+-ttf5flf)si"[2*2H=ose2] 


ri   r,  cos 8 o  ^  riJ  r,  cos8n    r~,  „    «  -, 

[1-  — L  rAj-[l+  — L  r^  cos[2k,Hcos8i] 

^   r   cose  ;  *»    r  cose  J  z     z 


(2.29) 


,3  (1  -^H|ff)sin(2X2Hcose2) 
-  tan  ±   ■ —J 

C1  +  r1  ^se^COS^2k2Hcos92) 
2       1 


In  the   experiment  to  be  described,   the  following  para- 
meters were  used  in  Equation  2.29: 

ri  =  p,c,  =  [998  £f]  [1481  m/s]  =  1.48  x  106  ^l3- 

x     ■  l  mJ  m 

r2  =  p2c2  =  [1200  £2]  [1992  m/s]  =  2.39  x  106  £2l£ 

*  L    ■  mJ  m 


H  »  3.18  x  10"3 


m 
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From  the  geometry  of  the  setup, 


8x  =  tan"1 


'rcvr 


^rcvr  +  dsourc 


=  tan_1  7*!*  -  «.3« 


where  <3rcvr  and  dsource  are  the  receiver  and  source 
depths,  respectively.  A  plot  of  the  phase  angle  caused 
by  the  plate,  Equation  2.29,  as  a  function  of  frequency 
with  the  parameters  shown  in  Figure  2.4  above  is  in 
Figure  2.4.  Note  the  conventional  of  a  negative  angle 
producing  a  phase  delay  in  the  signal. 

In  the  frequency  range  of  interest,  f  =  0  to  100 
kHz,  t  decreases  from  0°  at  f  -  0Hz,  to  -104*  at 
f  ■   100,000  kHz. 
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Figure   2.4      Phase   Delay  Caused   By   the   Plate 
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The  phase  delay  caused  by  the  plate  will  be  mani- 
fested by  a  shift  in  the  frequency  of  the  nodes  and 
antinodes  of  the  interference  pattern.  It  can  be  shown 
that  the  frequency  shift  Af  caused  by  the  phase  delay  x 
is : 

*£--£&  (2-30) 

where 

hr   =   difference  in  path  lengths  between  the  direct 

and  reflected  signals 
t   ■   phase  delay  caused  by  the  plate,   as  obtained 

from  Equation  2.29, 
c0  ■   speed  of  sound  in  the  fluid. 

G.   DESCRIPTION  OF  EXPERIMENT 

Two  small  LC-10  transducers  were  placed  in  the 
water,  one  acting  as  a  source,  the  other  as  a  receiver, 
and  were  connected  to  a  Hewlett  Packard  4162A  frequency 
analyzer  as  shown  in  Figure  2.5.  A  periodic  chirp  sig- 
nal of  0-100  kHz  was  applied  to  the  source  transducer, 
and  the  received  response  was  measured  and  plotted.  The 
LC-10  transducers  were  assumed  to  have  a  flat  response 
throughout  the  frequency  range.  While  the  response 
curves  of  the  LC-10  show  the  response  to  be  fairly  flat 
up  to  100-kHz,  the  assumption  is  really  not  tested, 
since   the    response   of   the    transducers   themselves 
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is  not  being  tested,  rather  the  difference   in  response 
with  and  without  the  pressure  of  the  plate. 

The  response  of  the  receiving  transducer  showed  the 
maxima  and  minima  in  the  pressure  at  its  location  due  to 
the  boundary.  The  response  was  measured  for  the  water- 
air  boundary,  and  without  changing  the  location  of  the 
two  transducers,  the  plate  was  introduced  and  the 
response  was  remeasured. 

H.   EXPERIMENTAL  RESULTS 

The  results  of  the  experiment  can  be  seen  from 
Figures  2.6  and  2.7.  Figure  2.6  is  a  plot  of  the  fre- 
quency response  of  the  receiving  hydrophone  with  fixed 
geometry  with  no  plate  present;  the  classic  Lloyd's 
mirror  experiment.  Note  the  frequency  pattern  of  the 
nodes  and  antinodes.  When  the  plate  was  placed  in  posi- 
tion without  moving  source  or  receiver,  the  positions  of 
these  nodes  and  antinodes  shift  to  the  right,  higher  in 
frequency,  as  can  be  seen  in  Figure  2.7.  This  was  pre- 
dicted as  the  manifestation  of  the  phase  angle  caused  by 
the  plate.  The  change  in  the  frequency  of  the  nodes 
increases  linearly  with  frequency,  corresponding  to  the 
increased  phase  angle. 
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Figure  2.5.   Lloyd's  Mirror  Experiment  set-up 
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Table  2.1  shows  a  comparison  of  the  predicted  nodal 
shift  to  the  nodal  shift  obtained  experimentally.  t  was 
obtained  from  Figure  2.4.  The  predicted  frequency  shift 
was  obtained  from  Equation  2.30. 


TABLE   2 . 

.1 

^node 

X 

-42° 

A*predicted 
1  kHz 

Afexp. 
1.02  kHz  ± 

37  kHz 

.13 

56 

-62 

1.5 

2.01      ± 

.43 

75 

-81 

1.9 

2.24      ± 

.37 

93 

-100 

2.3 

2.50      ± 

.39 

The  data  of  Table  2.1  are  shown  in  Figure  2.8.  The 
size  of  the  error  bars  of  the  experimental  data  are  due  to 
the  inaccuracies  in  reading  the  frequencies  of  the  nodes . 
Note,  however,  that  in  all  cases  the  experimental  data 
agreed  with  prerdicted  values  of  the  frequency  shift  with- 
in the  accuracies  of  measurement,  so  on  a  qualitative 
level,  at  least,  the  theoretical  predictions  of  the  nodal 
shifts  are  supported. 

In  further  comparison  of  Figures  2.6  and  2.7,  it  can 
be  seen  that  although  the  pattern  has  been  shifted  in 
frequency,  there  is  no  change  in  the  amplitudes  of  the 
maxima  and  only  a  small  change  in  some  of  the  minima  of 
the   frequency    response   pattern    when   the    plate   is 
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introduced.  This  suggests  that  there  is  insignificant 
change  in  the  amplitude  of  the  reflected  signal,  as  pre- 
dicted by  the  reflection  coefficient  R  having  a  magnitude 
of  one.  From  this,  it  is  inferred  that  there  is  insignif- 
icant translation  of  the  incident  wave  into  plate  wave 
modes  other  than  longitudinal .  The  ray-method  used  in  de- 
riving the  reflection  coefficient  R  assumes  that  the  com- 
pressional  wave  incident  at  the  plate  boundary  remains  a 
comressional  (or  longitudinal)  wave  in  the  plate. 

A  change  in  the  mode  of  the  wave  from  compress ional 
to,  for  instance,  flexural  or  shear,  would  involve  a  loss 
of  energy  in  the  reflected  wave,  which  would  be  seen  as  a 
change  in  the  amplitude  of  the  maxima  and  minima  of  the 
frequency  response  pattern.  Since  very  little  change  was 
observed,  it  appears  that  the  incident  wave  remained  a 
compressional  wave  as  it  traversed  the  plate  and  re- 
entered the  fluid  below. 

Thus,  it  has  been  demonstrated  that  the  effect  of 
placing  a  plate  into  a  fixed-geometry  method-of-images 
problem  is  to  introduce  only  a  phase  delay  into  the  re- 
flected signal,  causing  a  related  upward  shift  in  the  fre- 
quencies of  maxima  and  minima  in  the  pressure  for  the 
fixed  source  and  receiver  positions. 
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FREQUENCY  SHIFT  DUE  TO  PLATE 
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Figure   2.8      Frequency   Shift   due   to  Plate 
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III.   RADIATION  FROM  A  POINT  DRIVEN 
FLUID-LOADED  PLATE 

A.  INTRODUCTION 

Chapter  II  treated  the  effects  on  the  sound  field 
from  a  source  in  the  water  when  a  plexiglass  plate  is 
placed  on  the  surface.  This  chapter  will  deal  with  the 
acoustic  radiation  from  a  plate  driven  by  a  source 
attached  to  it.  When  the  source  is  driving  the  plate,  it 
excites  the  "acoustic  mode"  of  sound  propagation  in  the 
plate,  as  described  in  Chapter  I.  From  the  previous 
discussion,  the  acoustic  mode  leads  to  radiation  away 
from  the  plate  caused  by  the  spatial  gradient  of  the 
amplitude  of  the  flexural  wave  excited  by  the  source. 

The  principal  references  for  the  theoretical  de- 
scription of  the  problem  are  Feit  [Ref.  7]  and  Ross 
[Ref.  21].  A  computer  program  was  written  to  compute 
the  radiation  pattern  from  the  driven  plate  at  various 
frequencies .  An  experiment  was  conducted  to  obtain 
actual  radiation  patterns  from  a  plexiglass  plate  placed 
on  the  water. 

B.  THEORETICAL  DEVELOPMENT 

As  was  seen  in  the  background  section  of  Chapter 
One,  the  theory  concerning  the  point-driven  fluid-loaded 
plate  is  well  known,   and  its  full  development   will  not 
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be  given  here.  However,  for  continuity,  key  points 
leading  to  an  understanding  of  the  solution  will  be 
summarized. 

1.   Radiation  from  a  Propagating  Flexural  Wave 

A  flexural  wave  propagates   in  a  plate  with  wave 
number  kf  where 


1   cb 


(3.1) 


When  this  plate  is  in  contact  with  a  fluid,  vibrations 
of  the  plate  will  be  imparted  to  the  fluid,  producing 
waves  in  the  fluid  of  wave  number  k0  where 


ko  =  kl  +  kl 


(3.2) 


where  k^  and  k|   are  the  wave  numbers  perpendicular  and 
v parallel  to  the  plane  of  the  plate,  respectively. 


Figure  3.1   Sound  Radiation  from  a  Flexural 

Wave  of  a  Plate 
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The  length  of  this  vector  is  given  by 
k02  =  k±2  +  k„2  (3.3) 

the  component  parallel  to  the  plate  is 


V 


k02  -  k±2  (3.4) 


and  from  Snell's  Law  this  must  be  equal  to  kf  •  This 
leads  to  [Ref.  20,  p.  163] 

k±  =  JkQ2  -  kf2  (3.5) 

From  the  boundary  conditions  at  the  interface,  the 
normal   components  of  particle   velocity  must   be  equal, 

v,   =  u.  (3.6) 

Lf  io 

where   v^    refers  to  the  perpendicular  component  of  the 

plate  particle   velocity  and   uL   is   the  fluid  particle 

o 

velocity  perpendicular  to  the  plate. 

From  these  conditions,  it  can  be  shown  that  [Ref.  21 
p.  163]  the  acoustic  pressure  £  on  the  plate  due  to 
the  fluid  is: 

£  =  -jo)Poc0     kn    =  w  (3.7) 

^o2  "  kf2 

where  w  is  the  normal  displacement  of  the  flexural 
wave.  Equation  3.7  gives  the  acoustic  pressure  £  on 
the  plate  attributable  to  fluid  loading.  Note  that  £  is 
dependent  on  the  relative  values  of  k0  and  kf  . 
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This  leads  to  the  important  effect  of  coincidence. 
When  kf>k0  »  t*ie  expression  for  the  fluid  pressure  is 
real,  which  represents  mass-loading  [Ref.  21  p.  163]  and 
acoustic  power  is  not  radiated  into  the  fluid.  When 
k  =  kf  ,  the  condition  of  coincidence  is  reached.  The 
expression  for  the  radiated  acoustic  pressure  from  the 
flexural  wave  would  seemingly  predict  infinite  pressure. 
(This  cannot,  however,  occur,  because  fluid  loading  pre- 
vents the  two  speeds  from  being  exactly  equal  [Ref.  21 
p.  163])  When  kf<kQ  ,  the  expression  for  the  fluid  pres- 
sure in  Equation  3.7  is  imaginary,  so  the  effect  is  that 
of  a  radiation  resistance,  [Ref.  21  p.  163]  and  the  pro- 
pagating wave  radiates  power  away  from  the  plate. 

The  coincidence  frequency  is  obtained  by  equating 
the  wave  numbers  of  the  flexural  wave  in  the  plate  and 
the  compressional  wave  in  the  fluid  and  solving  for  fre- 
quency wc  .  But  because,  the  flexural  wave  speed  varies 
with  frequency  and  formulae  for  the  flexural  wave  speed 
depend  on  whether  thin  plate  or  thick  plate  approxima- 
tions are  used,  there  are  different  forms  for  the  solu- 
tion for  the  coincidence  frequency. 

As  there  were  three  approximations  for  the  flexural 
wave  speed  (Equations  1.6  -  1.8),  there  are  three 
formulae   for  the   coincidence   frequency.    In   the  low 


55 


frequency  limit,  or  thin-plate  theory,  the  coincidence 
frequency  is  given  by: 

«ci  =   X\   c°2  (3*8a) 

In  the  intermediate  frequency  range,  the  coincidence 
frequency  cannot  be  easily  solved  for  implicitly,  but  is 
given  by  the  solution  of: 


,     "ci.,.1/2    uci 

~i     I'  *  (— )')      -  '-5-  .  a!  U.8B) 


In   the  high   frequency  range,  the  coincidence  frequency 
can  be  given  by: 


c, 


12        c  2 


*ch  *  ^r^       2     S — 2  (3-8c) 

CP   "  2   °° 

Note  here   that   if   cp  <  Jj   cQ  ,    Equation   3.8c    is 

negative,  which  implies  that  coincidence  cannot  be 
reached.  Recall  that  the  high  frequency  asymptote  of 
the  flexural  wave  speed  is  the  Rayleigh  wave  speed,  a 
value  somewhat  less  than  the  shear  wave  speed.  If  this 
speed  is  less  than  the  speed  of  sound  in  the  fluid,  the 
condition  of  k0  =  kf,  or  coincidence,  cannot  be  reached. 
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2 .   Radiation  from  a  Point.  Source 

Heretofore,  we  have  discussed  the  flexural  wave 
itself  radiating  power  from  the  plate.  The  conclusion 
that  no  power  is  radiated  below  coincidence  applies  only 
to  infinite  plane  bending  waves  on  an  infinite  plate. 
If  excited  by  a  localized  force  even  below  coincidence, 
sound  will  be  radiated  because  the  waves  in  the  plate 
have  divergence  [Ref.  20  p.  169]. 

The  pressure  at  a  field  point  (R,  $ ,  9)  from  a 
point-driven  fluid- loaded  plate  can  be,  in  a  derivation 
by  Feit  [Ref.  7,  p.  1492],  eventually  expressed  in  the 
form  (3.9) 

£(R,*,8)-  ^^(H^^koRsin^sinwJe^o^^^^^vCwJsinwdw 

where 

p(R,<j),8)  =  acoustic  pressure  in  spherical 

coordinates  with  the  source  at  the 
origin, 

R  =  the  slant  distance  from  the  source, 

4>  =  the  azimuthal  angle, 

9  =  the  circumferential  angle  measured  in 
the  plane  of  the  plate, 

w  =  the  variable  of  integration, 

k0  =  the  wave  number  in  the  fluid, 

C  =  contour  of  integration  shown  in 
Figure  3.2, 

FQ  =  the  amplitude  of  the  source. 
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(1) 

H     =  the   zero-order   Hankel   function   of  the 
°      first  kind,  and 

V(w)  =  a  complicated  formula  that  relates  the 
vibrating  field  of  the  plate  to  the  water 
through  the  Mindlin-Timoshenko  thick 
plate  equation.  V(w)  contains  the  poles 
of  the  integral  in  Equation  3.10. 

The  solution  of  the  integral  was  obtained  by 
Brekhovskikh  [Ref.  18,  pp.  242-258]  and  is  given  by  Feit 
[Ref.  7,  Equation  24].  Ross,  however,  gives  a  more  use- 
ful form  of  Feit ' s  solution: 


p(R.»)-  VS^^   X  <3-10> 


_U)  2  Co  2  Cq  2 

cos4>    U„l     fsin2(j>-fcr]     Usin2<t)-rr  f  cp]    J 

"^       -*■*  ^Tc7r7^T^^~T~"  re,  —  * 


-l 


where  rp  is  the  plate  shear  parameter,  given  by 

rp  -  2.65(1+  |  a   +  |  a2)  (3.11) 

and   0   is   the  fluid   loading  factor,  given  by  Equation 
3.13,  which  is  discussed  fully  in  the  next  section. 

The  bracketed  term  in  Equation  3.10  along  with 
the  cos<f>  can  be  described  as  a  directivity  function 
H  (o>,<fr).  This  directivity  function  was  encoded  in 
FORTRAN  and  the  program  (given  in  Appendix  A)  was  used 
to  generate  plots  of   theoretical   directivity   patterns 
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at  several  different  frequencies.  From  these  figures 
(Figures  3.5-3.11)  the  coincidence  effect  can  be  readily 
observed.  The  figures  were  based  on  thin  plate  theory, 
and  cover  frequencies  from  below  to  above  coincidence. 
They  represent  a  hypothetical  plate  and  are  presented 
for  a  discussion  on  the  effects  of  coincidence  and  fluid 
loading.  They  will  be  discussed  in  section  El  of  this 
chapter. 
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Figure  3.2   Integration  Path  of  Integral  of  Equation  3.9 

(from  Feit) 


C.   THE  EFFECT  OF  FLUID  LOADING 

Although  Equation  3.10  is  valid  for  all  ranges  of 
plate  and  fluid  parameters  and  frequencies,  the  expres- 
sion may  be  simplified  greatly  and  further  insight  may 
be  obtained  by  considering  the  effect  of  fluid  loading. 
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The  fluid  loading  factor,   3  ,   is  the  ratio   of  the 

specific  acoustic   impedence  of   the  fluid   to   the  mass 

reactance  per   unit  are  of  the  plate.    Ross  [Ref.  20  p. 

166]  defines  it  as 

Po  co  / 0  , 0  n 

3  =  r—  (3.12) 

u>Ppt 

where 

pQ  =  density  of  the  fluid 

c0  -  speed  of  sound  in  the  fluid 

a)   =  radian  frequency  of  the  flexural  wave 

pp  =  density  of  the  plate 

t   =  thickness  of  the  plate 
Applying  the  formula  for  the  thin  plate  coincidence 
frequency,   Equation   3.8a,   the   fluid   loading   factor 
becomes 

1   P0  Cp  a»c 
3   ■  -^=  —  7T-   —  (3.13  ) 

vn  pp  co  w 

where 

cp  ■  the  plate  longitudinal  wave  speed,  given  in 
Figure  1.1a, 

u)c  ■  the  coincidence  frequency. 

A  fluid  loading  factor  much  less  than  1  corresponds 
to  light  fluid  loading,  which  occurs,  for  example,  in  a 
metal  plate  in  air.  A  fluid  loading  factor  much  greater 
than  1  corresponds  to  heavy  fluid  loading  found,  for 
example,  a  plexiglass   plate  on  water.    The   effect   of 
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fluid  loading  can  best  be  seen  by  considering  four 
general  cases:  frequency  below  coincidence,  light  and 
heavy  fluid  loading;  frequency  above  coincidence/  light 
and  heavy  loading. 

When  u)<o)c  and  3<<1,  a  case  of  light  fluid  loading, 
the  radiation  pattern  is  that  of  a  small  baffled  piston, 
or  monopole-type  radiation,  non-directive.  For  this 
case,  Equation  3.10  reduces  to  [Ref.  7  p.  1492] 

P(r.t.e)  =£r°  ejkoR  (3.14) 

Ross  [Ref.  20  p.  170]  states  that  the  radius  of  this 
piston  is  the  sum  of  the  radius  of  the  surface  over 
which  the  force  is  applied  and  a  quarter  of  the  flexural 
wavelength.  Ver  and  Holmer  [Ref.  12  p.  278]  give  a 
slightly  different  form.  They  state  that  for  a  plate 
driven  by  point  source,  the  equivalent  radius  of  the 
baffled  piston  is  0.286  \c  ,  where  Xc  is  the  flexural 
wavelength  at  coincidence. 

When,  o)<o)c  and  3>>1  ,  heavy  fluid  loading,  the  point 
force  acting  on  a  thin  plate  radiates  sound  as  through 
the  plate  were  not  there,  representing  dipole-type 
radiation,  and  the  radiation  pattern  is  cpsinusoidal . 
Equation  3.10  reduces  to  [Ref.  7  p.  1492]: 

P(R'*)  =  "2TIT2  eJk°R  cos+  (3'15) 
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When  o)>o)c  and  &<<1  ,  the  radiation  pattern  is   still  co- 
sinusoidal  [Ref.  21  p.  171]. 

When  u)>wc  and  3>>1  ,  Equation  3.10  can  be  reduced 
only  slightly,  and  gives  the  interesting  radiation  pat- 
terns corresponding  to  the  flexural  wave  radiation  and 
radiation  from  the  driving  piston. 

D.   THE  EFFECT  OF  STRUCTURAL  ABSORPTION1* 

Above  coincidence,  then,  the  absorption  of  the  flex- 
ural wave  as  it  propagates  through  the  plate  will  effect 
the  radiation  pattern.  Feit  [Ref.  7,  p.  1494]  states 
this  and  states  that  (without  proof)  the  far  field  pres- 
sure maximum  is  reduced  by  5  dB  for  a  loss  factor, 
ri=0.1.  (See  Equation  3.18  for  the  definition  of  the  loss 
factor  n  .)  Ungar  [Ref.  22,  p.  453]  lists  a  loss  factor 
typical  of  plexiglass  of  n  =  0.02  to  0.04  in  the  audio 
frequency  range.  Therefore,  it  is  expected  that  the  de- 
crease in  predicted  peak  pressure  amplitude  due  to  the 
absorption  of  the  flexural  wave  will  be  less  than  5  dB. 

An  experiment  was  conducted  to  quantify  the  loss 
factor  for  the  plate  used  in  the  experiment.  This  was 
done  also  to  investigate  the  effects  of  the  reflection 
of  the  flexural  wave  from  the  edge  of  the  plate,  i.e., 


**Although  the  effect  is  the  same,  the  term  absorption 
is  perhaps  more  correct  than  damping  as  the  term  damping 
implies  additional  materials  applied  to  the  plate. 
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patterns  were  induced.  (Other  samples,  for  instance 
of  metalj  did  produce  Schladni  patterns  using  this 
technique. ) 

Qualitatively,  the  effect  of  the  plexiglass  loss 
factor  n  was  very  apparent.  Refer  to  Figure  3.3  which 
is  a  plot  of  flexural  wave  attenuation,  after  correction 
for  cylindrical  spreading.  As  the  accelerometer  was 
moved  farther  away  from  the  source,  the  corrected  ampli- 
tude of  the  flexural  wave  steadily  decreased  with  range, 
dropping  40  dB  by  the  end  of  the  plate  for  the  frequen- 
cies of  interest.  No  reflections  of  the  pulse  off  the 
edges  were  observed. 

Ver  and  Homer   [Ref.  12,  p.  276]   give  the  losses  in 
power  radiated   from  the   plate  due   to  absorption  to  be 
e-nkf  R/2 

where     kf  =  flexural  wave  number  (3.18) 

R   ■  range  from  rource 
n   =  loss  factor  of  the  plate  material 

using  this,  with  the  100  kHz  absorption  curve  (the  only 
one  straight  enough  to  be  trusted)  yields  n  ~  0.03 
which  is  consistent  with  the  values  of  Ungar  [Ref.  21, 
p.  453]. 

According  to  Junger  and  Feit  [Ref.  13,  p.  249],  the 
effects  of  structural  damping  on  the  far-field  pressure 
pattern  are  accounted  for,   for  small  loss   factors,   by 
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multiplying  c_  in  Equation  -3.10  by  (l-jn)1/2 
-  (1-1/2  jn)#  or  wc  by  (1+  1/2  jn).  They  further  state 
that  damping  plays  no  significant  effect  on  the  maximum 
pressure  of  the  radiation  field  below  coincidence.  This 
is  also  intuitively  true  because  below  coincidence  the 
propagating  flexural  wave  does  not  radiate  and  therefore 
its  absorption  has  no  effect.  Above  coincidence,  they 
state  (as  Feit  did  in  his  earlier  work)  that  the  effect 
is  to  reduce  the  far-field  pressure  maximum.  Figure 
3.4,  taken  from  Feit  [Ref.  7,  p.  1494],  show  the  effect 
of  structural  damping  on  the  far-field  of  a  point- 
excited  thin  plate. 


40°  50" 


60°  70°  80° 

Poiar    Angle   ♦ 


Figure  3.4,   The  Effect  of  Damping  on  the 

Far-Field  Pressure  (from  Feit) 
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E.   DISCUSSION  OF  THEORETICAL  PLOTS 

1 .   The  General  Effects  of  Increasing  Frequency 

Figures  3.5  to  3.11  were  generated  from  the  pro- 
gram given  in  Appendix  A.  They  show  the  radiation 
directivity  pattern  H( $) ,  normalized  to  the  on-axis 
pressure,  for  a  frequency  range  from  below  coincidence 
to  above  coincidence. 

The  plate  parameters  used  in  these  figures  were 
those  of  a  hypothetical  plate  and  fluid  combination. 
They  do,  however,  show  quite  well  the  effects  of  in- 
creasing frequency  on  the  radiation  pattern  in  the 
fluid.  Thin-plate  theory  was  assumed,  and  the  ratio  of 
oj/o)c  was  the  parameter  that  was  varied. 

Figures  3.5  to  3.7  show  the  directivity  patterns 
below  coincidence.  Far  below  coincidence, Figure  3. 5, the 
radiation  from  the  plate  is  cosinusoidal,  with  maximum 
radiation  at  |  =  0°  and  minimum  radiation  at  $  =  90°. 
It  was  shown  in  Section  Bl  of  this  chapter  that  below 
coincidence  an  infinite  plane  propagating  flexural  wave 
will  not  radiate  sound.  It  was  shown  in  Section  C  that 
the  radiation  pattern  should  be  that  of  a  low  frequency 
small  radius  piston  mounted  on  an  infinite  baffle.  This 
is  revealed  in  these  first  three  figures.  Note  that  as 
the  coincidence  frequency  is  approached,  the  pattern 
changes  into  a  wider  (Figures  3.7  and  3.8)  pattern. 
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As  the  frequency  of  the  source  becomes  larger 
than  the  coincidence  frequency,  the  phase  speed  of  the 
propagating  flexural  wave  exceeds  that  of  the  fluid,  and 
therefore  the  propagating  flexural  wave  can  radiate 
energy  away  from  the  plate.  This  is  seen  in  Figures  3.9 
and  3.10. 

Well  above  coincidence,  Figure  3.11,  a  third 
lobe  is  predicted.  This,  according  to  Feit,  is  the 
"...  nonpropagating  near-field  deflection  of  the 
plate  ....  itself  becoming  a  propagating  wave." 
[Ref.  7,  p.  1493] 

Feit,  in  a  later  work  with  Junger  [Ref.  13,  p. 
248],  state  that  thick  plate  theory  yields  larger  values 
for  the  angles  at  which  the  propagating  flexural  wave 
radiates  and  smaller  values  of  peak  pressure.  That  is, 
the  lobes  in  the  pattern  will  be  at  a  greater  angle  and 
have  a  smaller  value.  Figure  3.11,  from  Feit  and 
Junger,  show  a  comparison  of  the  size  of  the  pressure 
lobe  and  the  angle  at  which  it  appears  for  thin-plate 
theory  and  thick-plate  theory.  The  figure  neglects  the 
effect  of  structural  damping. 
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Fig.  3.11.  Comparison  of  Size  and  Angle  of  Pressure 
Lobes,  Thin-Plate  Theory  vs.  Thick-Plate 
Theory  (From  Junger  and  Feit) 


2.   The  Effect  of  the  Fluid  Loading  Factor,  B 

The  plots  of  the  previous  section  showed  the 
general  trend  in  the  radiation  pattern  when  frequency  is 
increased  for  a  given  set  of  thin-plate  parameters. 
However,  there  are  subtle  changes  that  take  place  in  the 
patterns  when  different  parameters  are  used.  The  den- 
sity of  the  plate  can  be  measured  with  reasonable 
accuracy,  and  Poisson's  ratio  should  be  somewhat  similar 
among  plexiglass  samples.  However,  Young's  modulus  and 
the  exact  effect  of  Poisson's  ratio  in  the  wave  speed 
formula  may   vary.    Additionally,   Equation   3.10  lends 
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itself  better  to  changing  certain  parameters  than 
others .  The  parameter  chosen  for  variation  was  the 
fluid  loading  factor,  3  ,  for  fixed  ratios  of  u)C/o>  • 
Figures  3.13  to  3.21  are  the  results  of  these  runs.  The 
effect  of  varying  the  fluid  loading  factor  3  as  follows. 
Increasing  the  fluid  loading  smooths  the  radiation 
pattern,  reducing  the  peaks  and  making  the  coincidence 
effect  less  dramatic. 

Below  coincidence/  fluid  loading  should  not  have 
a  great  effect  on  the  radiation  pattern.  This  is  seen 
in  Figures  3.13  thru  3.15.  A  change  in  the  fluid  load- 
ing factor  3  for  a  frequency  ratio  of  <uc/u)  =  2  below  co- 
incidence has  a  minimal  effect  on  the  radiation 
pattern. 

At  coincidence,  the  effects  of  fluid  loading  on 
the  radiation  pattern  can  be  seen:  for  light  fluid 
loading,  Figure  3.16  shows  a  well-developed  radiation 
pattern,  but  as  3  is  increased  (Figures  3.17  and  3.18) 
the  radiation  pattern  is  less  pronounced  and  the  off- 
axis  pressure  gain  not  as  great. 

Similarly,  above  coincidence,  (coc/w)  ■  0.86,  as 
seen  in  Figures  3.19  through  3.21,  the  increase  in  fluid 
loading  causes  a  less  dramatic  radiation  pattern. 

The  importance  of  3  is  noticeable.  Even  though 
in  the  computer  runs   3  was   never  "much  less   than  one" 
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which  was  described  as  the  condition  for  light  fluid 
loading,  small  changes  in  3  caused  large  changes  in  the 
radiation  patterns . 

F.   DESCRIPTION  OF  EXPERIMENT 

An  experiment  was  conducted  to  measure  the  radiation 
pattern  from  a  point-driven  sheet  of  plexiglass .  See 
Figure  3.22.  A  l/8th-inch  thick  sheet  of  plexiglass  of 
dimensions  87cm  x  152cm  was  floated  on  the  surface  of  a 
freshwater  anechoic  tank.  A  small  silicon  bead  was 
applied  to  the  edges  of  the  plate  to  keep  it  afloat.  In 
the  center  of  the  plate  a  Tonpilz  longitudinal  vibrator 
of  the  author's  design5  and  manufacture  was  bolted  to 
the  plate.  The  vibrator  was  driven  by  an  HP-3314A 
function  generator  operating  in  the  pulse  mode.  A  maxi- 
mum driving  output  voltage  of  30V  was  applied. 

A  schematic  cross-section  of  the  Tonpilz  vibrator  is 
given  in  Appendix  B.  Its  tail  mass  was  made  out  of 
stainless  steel  and  its  head  mass  out  of  aluminum.  This 
was  to  impart  as  large  a  signal  as  possible  to  the  plex- 
iglass sheet  while  keeping  the  vibrator's  natural 
resonance  frequency  fairly  high.  An  admittance  circle 
of  susceptance  vs.  conductance  and  plots  of  admittance 


5For   an   excellent   description   of   the   design  of 
Tonpilz  vibrators,  see  Wilson  [Ref.  23,  Ch.  6]. 
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and  susceptance  vs.  frequency  for  the  vibrator  are  also 
given  in  Appendix  B.  They  show  the  motional  resonance 
frequency  of  the  vibrator  to  be  92.4  kHz,  and  a  Q  of 
about  112. 

It  is  stated  that  although  the  motional  resonance 
frequency  of  the  vibrator  was  within  the  frequency  range 
of  interest  for  this  experiment,  this  was  not  of 
key  importance  because  a  swept- frequency  mode  was  not 
used.  However,  what  is  of  interest  in  the  pulse  mode  is 
the  Q  of  the  vibrator.  The  Q  of  the  vibrator  showed  its 
effect  in  the  transient  response  of  the  vibrator  in  the 
pulse  mode.  The  transient  response  decay  time  grew  when 
the  driving  frequency  approached  the  motional  resonance 
frequency  of  the  vibrator.  However,  in  every  data  set 
measurement,  the  pulse  length  was  adjusted  to  give  maxi- 
mum steady  state  response  time  for  the  given  pulse 
interval . 

The  head  mass  of  the  vibrator  was  about  1  cm  in 
radius  and  was  attached  to  the  plate  with  a  thin  film  of 
silicon  grease  to  ensure  positive  contact  between  the 
head  mass  and  plate.  The  clamp  rod  of  the  vibrator  was 
used  to  compress  the  vibrator  against  the  plate. 

The  radiation  pattern  from  the  driven  plate  was  mea- 
sured with  a  Celesco  LC-10  hydrophone  whose  output  was 
sent  through  a  Tektronics  AM-502  low  noise  preamplifier 
to  a  Kikusui  oscilloscope.    The  LC-10  was  mounted  on  an 
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Figure  3.22,   Driven  Plate  Experimental  set-up 
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angular  jig  that  allowed  it  to  be  swept  from  <j>  =  -90°  to 
90°  at  constant  range  from  the  vibrator.  An  Endevco 
K920  accelerometer  was  used  at  times  to  probe  the  plate 
to  obtain  attenuation  and  edge-reflection  data. 

The  pulse  mode  was  used  because  of  a  problem  of 
electrical  pickup  between  the  leads  to  the  vibrator  and 
the  LC-10 .  The  electrical  pickup  could  not  be  removed 
even  with  the  use  of  grounded  coaxial  cable  and  shielded 
connectors . 

The  length  of  the  pulse  varied  with  frequency.  The 
limiting  factors  for  the  pulse  length  were  the  pulse  in- 
terval and  the  transient  response  time  of  the  vibrator. 
The  pulse  interval  was  limited  because  the  LC-10  was 
located  at  a  fixed  distance  from  the  source  and  recep- 
tion of  a  pulse  by  the  hydrophone  had  to  be  ended  before 
the  next  began  to  prevent  overlapping  signals.  The 
transient  response  time  of  the  vibrator  also  limited  the 
length  of  pulse  used.  The  received  pulse  showed  the 
transient  and  stead-state  response  of  the  vibrator. 
However,  amplitude  data  could  only  be  taken  on  the 
steady-state  portion  of  the  received  signal.  This 
further  limited  the  time  available  in  the  pulse  with 
which  to  drive  the  plate  in  a  steady-state  condition. 

Both  the  pulse  interval  and  pulse  length  were  ad- 
justed  for  each   data  set   measurement  to   provide   the 
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maximum  time  available  with  which  to  drive  the  plate  in 
a  steady-state  condition. 

The  properties  of  the  plexiglass  used  were  deter- 
mined experimentally  by  Johnson  and  Denny  [Ref.  23.] 
Table  3.1  lists  the  wave  speeds  for  the  three  wave  types 
they  measured.    Note  that  both  the  shear  wave  speed  and 


TABLE   3 . 1 

Experimental  Wave  Speeds 
for  Plexiglass 

CP  = 

2240  m/s 

cs  = 

1373  m/s 

cb  = 

700-1150  m/s 

therefore  the  high-frequency  asymptote  of  the  flexural 
wave  speed  are  less  than  the  speed  of  sound  of  the 
water.   Therefore,  coincidence  cannot  be  reached. 

G.   EXPERIMENTAL  RESULTS 

1 .   Discussion  of  Plots 

Several  plots  were  obtained  at  frequencies  from 
50  kHz  to  245  kHz  with  varying  degrees  of  success.  As 
was  stated  earlier,  the  high  frequency  asymptote  of  the 
flexural  wave  speed  was  less  than  the  speed  of  sound  in 
the  fluid,  so  the  coincidence  effect  could  not  be 
obtained.  Still  there  was  a  change  in  the  radiation 
pattern  from  the  plate  as  frequency  was  increased. 
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Figure  3.23  shows  the  experimental  radiation 
pattern  when  the  plase  was  driven  at  f=50  kHz.  There 
are  several  aberrations  in  the  pattern  when  the  deflec- 
tion angle  increased  above  40°  .  Two  plots  were  made, 
one  with  the  hydrophone  at  a  range  of  0 . 5  m  and  the 
other  at  a  range  of  1  m  .  This  was  done  to  see  if  the 
distance  from  the  source  to  the  vibrator  had  an  effect 
on  the  pattern.  It  can  be  seen  that  it  did  not,  as  the 
two  curves  appear  to  have  quite  similar  behavior. 
Still,  the  general  cosinusoidal  pattern  of  a  point- 
driven  plate  will  below  coincidence  is  seen. 

In  Figure  3.24,  the  aberrations  of  the  previous 
pattern  are  more  pronounced.  The.  next  sets  of  data, 
Figures  3.25  and  3.26  at  92.3  and  93.5  kHz,  respectively 
show  the  development  of  a  side  lobe  at  an  angle  of  about 
45°  .  Figure  3.26,  where  data  were  taken  at  angles  of 
-90°  to  90°,  show  an  asymmetry  in  the  pattern.  This  was 
seen  in  a  few  plots  and  is  thought  to  have  been  caused 
by  asymmetrical  excitation  of  the  plate.  This  is  dis- 
cussed further  in  section  G.3. 

The  next  three  plots,  Figures  3.27  to  3.29,  at 
frequencies  of  100  kHz,  120  kHz,  and  150  kHz  show  the 
development  of  lobes  at  about  40°.  There  are  multiple 
aberrations  in  the  pattern  which  may  or  may  not  have 
been  additional  minor  lobes.   Note  that  only  the  pattern 
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at  f  =  150  kHz  show  a  clean  major  and  minor  lobe 
pattern. 

In  the  following  three  sets  of  plots,  figures 
3.30  to  3.32  at  f  =  170  kHz,  194  kHz,  and  196  kHz,  the 
well-developed  pattern  of  Figure  3.29  is  lost  and  the 
pattern  became  more  complex,  with  multiple  aberrations 
which  may  have  represented  small  minor  lobes  in  the  pat- 
tern or  interference  resulting  from  incomplete  and  asym- 
metrical excitation  of  the  plate. 

The  final  three  data  sets,  figures  3.33  to  3.35 
at  frequencies  of  225  kHz,  227  kHz,  and  246  kHz,  respec- 
tively, show  patterns  that  are  very  complex,  and  as  evi- 
dence in  figures  3.33  and  3.35,  very  asymmetrical.  The 
side  lobes  of  the  pattern  are  just  as  large  as  the  on- 
axis  lobe. 

2 .   Error  Analysis 

The  method  of  data  collection  was  as  follows. 
For  a  given  frequency,  with  the  hydrophone  on-axis  (0°) 
the  output  voltage  of  the  preamplifier  was  adjusted  to 
7.0V.  All  readings  thereafter  were  immediately  corrected 
to  dB  re  7.0V.  The  signal-to- noise  ratio  (SNR)  varied 
between  data  runs  because  of  the  output  response  of  the 
vibrator.  On  the  average,  however,  the  SNR  was  esti- 
mated to  be  30  dB  at  zero  deflection  angle. 
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Since  the  data  were  normalized  to  the  on-axis 
reading,  the  error  increased  with  decreasing  sound  pres- 
sure (and  hence  voltage)  levels.  Table  3.2  shows  the 
resulting  dB  error  for  various  levels. 


EXPERIMENTAL  ERROR 

TABLE 

3.2 

SPL 

ERROR 

OdB  re  7V 

±  0.12 

dB 

-5 

±  0.56 

dB 

-10 

±0.84 

dB 

-15 

±  1.27 

dB 

Angular  measurements  were  quite  precise  with  the 
jig  that  was  used.  Angles  were  precise  to  ±0.1  degree. 
As  precise  as  the  angular  measurements  were,  they  were 
not  as  accurate.  While  every  attempt  was  made  to  accu- 
rately position  the  hydrophone  directly  underneath  the 
vibrator  at  the  start  of  each  data  collection  period, 
parallax  and  the  water  refractive  effect  hampered  accu- 
rate positioning.  It  is  estimated  that  the  angular 
accuracy  was  ±5  degrees .  Techniques  to  improve  the 
accuracy  were  abandoned  because  it  was  reasoned  that  the 
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the  result  of  any  inaccuracies  would  be  an  angular 
skewing  of  a  pattern,  not  a  change  in  the  pattern 
itself. 

It  would  seem  therefore,  that  inaccuracies  or 
errors  in  data  collection  were  not  the  cause  of  such 
large  aberrations  are  asymmetries  in  the  experimental 
plots. 

3 .   Discussion  of  Results 

Since  the  high  frequency  asymptote  of  the  flex- 
ural  wave  speed  was  less  than  that  of  the  speed  of  sound 
in  the  fluid,  coincidence  was  not  possible.  Therefore, 
the  large  off-axis  pressure  maxima  associated  with  the 
driven  plates  above  coincidence  were  not  obtained.  This 
may,  however,  be  proper  modelling  of  the  ice.  Stein 
[Ref.  2,  p.  288-289]  found  in  his  Arctic  ice  experiments 
that  " .  .  .  at  no  point  is  the  first  antisymmetric  mode 
[i.e.,  the  flexural  wave]  supersonic  with  respect  to  the 
water;  there  is  no  coincidence  frequency  above  which 
phase  speeds  are  faster  than  the  water  wave  speed.  This 
implies  that  in  the  plane-wave  limit  the  first  antisym- 
metric mode  never  radiates  and  never  loses  energy  to  the 
water . " 

After  a  careful  search  of  recent  literature  on 
driven  plates  below  coincidence,  no  theories  were  found 
that  would   suggest  a  radiation  pattern  other   than  that 
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of  a  cosinusoidal  pattern.  Nor  were  any  theories  found 
on  the  radiation  patterns  from  such  acoustically  slow 
plates  that  coincidence  was  not  possible.  Therefore, 
theoretical  plots,  such  as  those  developed  for  earlier 
sections,  were  not  available.  However,  some  possible 
explanations  are  offered  about  the  complexities  and 
asymmetries  of  the  patterns . 

Since  this  was  a  finite-sized  plate  it  could  be 
argued  that  plate  modes  of  vibration  were  set  up,  which 
caused  nodes  in  the  plate  that  would  affect  the  pressure 
field  in  the  fluid.  However,  the  pulse  mode  of  opera- 
tion was  used,  which  reduces  the  likelihood  that  this 
was  responsible  for  discrepancies.  Also,  it  was  found 
(and  stated  earlier  in  this  chapter)  that  no  reflections 
were  observed  from  the  plate  edge.  Therefore,  excita- 
tion of  normal  modes  of  vibration  in  the  plate  are 
thought  to  be  unimportant  for  the  pulse  experiment. 

Another  possible  explanation  is  that  the  hydro- 
phone was,  at  0.5  m,  in  the  "near- field"  of  the  plate. 
Feit  [Ref.  7,  p.  1482]  used  the  asymptotic  approximation 
kQR>>l  in  the  evaluation  of  the  integral,  Equation  3.9. 
At  the  lowest  frequency  at  which  data  were  taken,  50 
kHz,  kQR  equals  106,  so  it  would  appear  that  the  far- 
field  assumption  is  valid. 

The  complexities  of  the  patterns  might  also  be 
explained  by  the  pulse  not  being  long  enough.   The  pulse 
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interval  was  limited  to  0.30ms,  based  on  the  speed  of 
the  pulse  through  the  water  to  the  receiver  and  the 
transient  response  time  of  the  vibrator.  Table  3.3 
gives  values  of  travel  distance  that  a  flexural  wave 
packet  would  achieve  in  the  plate  in  0.30  ms  for  the 
frequencies  used. 


TABLE 

3 

.3 

f (kHz) 

c>^(m/s 

0 

d( 

m) 

50 

725 

0. 

22  ' 

100 

1026 

0. 

31 

200 

1150 

0. 

35 

The  cosinusoidal  patterns  expected  below  coinci- 
dence are  based  on  an  assumption  of  cylindrical  diver- 
gence of  the  flexural  wave  on  an  infinite  plate.  That 
is, the  pressure  at  a  field  point  in  the  fluid  is  made  up 
of  the  sum  of  pressures  caused  by  the  deflection  of  each 
incremental  area  of  the  plate.  If  the  entire  plate  were 
not  excited,  this  sum  would  be  incomplete.  From  Table 
3.2,  a  pulse  of  flexural  wave  of  frequency  50  kHz  will 
only  travel  out  0.22m  in  0.3  m/s.  Thus,  the  entire 
plate  was  not  being  excited.  Even  at  200  kHz,  the 
entire  plate  was  not  excited.  This  certainly  could  be 
one  of  the  contributing  factors  in  the  complexity  of 
experimental  plots. 
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Prior  to  the  discovery  of  the  electrical  pickup 
problem  mentioned  earlier  in  this  chapter,  the  plate  was 
driven  in  the  CW  mode.  The  complexities  in  the  patterns 
were  still  there.  This  would  indicate  that  use  of  the 
pulse  mode  was  only  a  contributing  factor,  not  the  main 
cause  of  the  complexities  of  the  experimental  plots. 

A  more  important  factor  is  thought  to  be  the  ab- 
sorption of  the  flexural  wave  in  the  plexiglass. 
Following  the  same  argument  above,  that  is,  the  pressure 
at  a  field  point  is  the  sum  of  the  pressures  caused  by 
the  individual  deflections  of  each  incremental  area  of 
the  plate,  the  effect  of  the  flexural  wave  absorption 
would  be  to  greatly  reduce  the  contributions  of  radia- 
tions from  the  incremental  areas  far  from  the  source . 

The  so-called  "acoustic  mode"  of  sound  propa- 
gation, which  this  experiment  was  modelling,  assumes 
sound  radiation  away  from  the  source  in  the  fluid  from  a 
lossless  plate.  Plexiglass  is  a  highly  absorptive 
medium;  therefore,  the  propagating  flexural  wave  decayed 
in  amplitude  due  to  losses  as  well  as  divergence  as  it 
travelled  away  from  the  vibrator.  Consequently,  the 
incremental  areas  of  the  plate  far  from  the  source  did 
not  deflect  nearly  as  much  as  those  close  to  the  course, 
so  they  did  not  contribute  as  much  to  the  pressure  at  a 
given  point  in  the  fluid. 
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The  plate  might  be  thought  of  as  a  complex  aper- 
ture, or  a  quasi-piston,  with  a  radius  of  the  head  mass 
radius  and  a  portion  of  the  adjacent  plate.  A  rough 
estimate  of  the  aperture  radius  is  the  head  mass 
radius  plus  either  one-fourth  the  flexural  wave  radius 
plus  either  one-fourth  the  flexural  wave  radius  or  one 
plate  thickness .  Both  give  similar  results .  Based  on 
this,  kQa  for  radiation  into  the  fluid  varies  from  2.79 
for  a  frequency  of  50  kHz  to  14  at  f  =  250  kHz. 

Even  though  the  aperture  is  shaded  and  phased, 
it  is  noteworthy  that  for  a  kQa  in  excess  of  4,  a  piston 
is  beginning  to  develop  a  side  lobe,  and  for  a  kQa  of 
10,  it  has  about  two  side  lobes.  This  is  not  inconsis- 
tent with  the  angular  scale  of  the  lobelike  structure 
observed. 

The  asymmetry  seen  in  a  few  of  the  radiation 
patterns  might  have  been  caused  by  radially  non-uniform 
excitation  of  the  plate.  As  was  mentioned  in  section  F 
of  this  chapter,  the  vibrator  was  compressed  against  the 
plate  by  a  bolt.  However,  asymmetrical  compression  of 
the  vibrator  against  the  plate  was  possible.  Possible 
causes  of  asymmetrical  compression  are  improper  machin- 
ing of  the  head  mass,  imperfect  centering  of  the  bolt 
through  the  vibrator  assembly,  or  a  non-uniform  layer  of 
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grease  between  the  head  mass  and  plate.  A  non-uniform 
radial  excitation  of  the  plate  seems  to  be  a  plausible 
explanation  for  asymmetry  in  the  radiation  patterns . 

Two  important  points  that  can  be  made  are  that 
it  appears  that  the  radiation  pattern  is  a  function  of 
the  source,  not  any  flexural  wave  propagation. 
Secondly,  situations  where  coincidence  is  not  possible, 
such  as  a  plexiglass  plate  on  water,  present  radiation 
patterns  at  lower  frequencies  similar  to  radiation  pat- 
terns below  coincidence  in  situations  where  coincidence 
is  possible. 
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IV.   CONCLUSIONS 

This  study  was  a  low-cost,  "low-tech",  first-cut 
attempt  to  duplicate  a  sound  source  mounted  on  the  ice 
in  the  Arctic.  The  first  experiment  performed  showed 
that  when  a  plate  is  introduced  at  the  boundary  of  a 
Lloyd's-mirror-type  problem,  the  effect  of  the  plate  was 
only  to  introduce  phase  delay,  causing  a  shift  in  the 
node  and  antinode  pattern.  If  any  flexural  waves  were 
excited  in  the  plate,  which  might  have  caused  an  in- 
crease in  amplitude  in  the  received  signal,  they  would 
have  been  absorbed  by  the  plexiglass. 

In  the  second  experiment  performed,  it  was  found 
that  even  though  the  theory  behind  a  point-driven  fluid- 
loaded  plate  is  well  known  and  many  interesting  results 
can  be  gleaned  from  that  theory,  extracting  experimental 
plots  upon  which   a  theory   can  be   based   is  difficult. 

Plexiglass  was  originally  chosen  because  its  speci- 
fic acoustic  impedance  is  fairly  close  to  both  ice  and 
water.  However,  the  flexural  wave  absorption  suffered 
in  commonly  available  plexiglass  and  more  importantly, 
its  low  flexural  wave  speed  render  it  a  poor  choice  to 
use  in  the  lab  in  this  type  of  experiment,  where  wave 
generation  and  propagation  are  important. 
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Several  interesting  types  of  plastics  are 
available.  Often  the  manufacturer  lists  the  density  and 
Young's  modulus  of  each  in  its  catalog.  However,  flex- 
ural  wave  attenuation  data  are  never  available.  It  is 
recommended  that  the  next  step  in  research  along  these 
lines  be  to  collect  samples  of  each  type  of  plastic  and 
measure  these  parameters.  The  time  available  for  this 
study  prevented  such  an  analysis. 
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APPENDIX   A 
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APPENDIX  B 
DATA  ON  TONPILZ  VIBRATOR 
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(35.5  g) 
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(6.0  g) 
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7/8"    > 


Figure  B.l.   Diagram  of  Tonpilz  vibrator 
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Figure  B.2.      Susceptance   and  Conductance  vs 

Frequency   of   Tonpilz   vibrator 
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